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Abstract — Optimal search strategies where targets are ob- 
served at several different angles are found. Targets are 
assumed to exhibit rectangular symmetry and have a uniformly- 
distributed orientation. By rectangular symmetry, it is meant 
that one side of a target is the mirror image of its opposite 
side. Finding an optimal solution is generally a hard problem. 
Fortunately, symmetry principles allow analytical and intuitive 
solutions to be found. One such optimal search strategy consists 
of choosing n angles evenly separated on the half-circle and 
leads to a lower bound of the probability of not detecting 
targets. As no prior knowledge of the target orientation is 
required, such search strategies are also robust, a desirable 
feature in search and detection missions. 

I. INTRODUCTION 

In mine hunting operations it is known that the detection 
performance improves when a target is observed many times 
at different aspect angles [l]-[6]. Similarly, classification 
algorithms [7]— [9] and sensor-arrays deployed for target 
localization and tracking [10]— [13] benefit from multi-aspect 
observations. This fact is, however, often overlooked. For 
example, the formula for the probability of detecting a target 
in a random search derived by Koopman is widely used yet 
it assumes no angular dependence [14]. 

In this paper, search strategies that minimize the overall 
probability of not detecting a target observed at several 
different angles are identified. Finding an optimal angular 
configuration is a priori intractable as it is multi-dimensional 
in the sense that each observation is independent of one 
another and hence each observation angle must be considered 
as a separate dimension. What is more, the explicit expres- 
sion for the overall probability of detection can be hopelessly 
complicated even when the probability of detection for a 
single observation is simple and the number of observations 
is small. 

The novelty of our approach lies in the fact that this 
normally intractable problem is solved using an elegant 
symmetry argument. Specifically, targets are assumed to 
exhibit rectangular symmetry. That is, the left-hand side of a 
target is the mirror image of its right-hand side, and its rear 
end is the mirror image of its front end. Many targets can be 
approximated with this class of symmetry including canoes, 
ships, submarines, mines and human bodies. 

Optimal angles are then shown to be evenly distributed 
on multiples of the half-circle as in Ref. [13]. However, this 
constitutes a departure from the current literature on sensor 
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Fig. 1. Rectangular target observed at angle x. 



geometry [10]— [13] as our result is derived for finite-extent 
targets rather than for point targets. The simplicity of the 
solution implies that no complicated calculations are required 
prior to a search as long as the target has the assumed 
approximate symmetry. This fact should improve the task 
of planning the path of mobile sensors, such as unmanned 
vehicles, to search for fixed targets, as well as of deploying 
a fixed sensor array to monitor traffic through choke points. 

Assumptions and the minimization problem are stated in 
Section|n] while Sectionlllllpresents a set of search strategies 
that minimizes the probability of no detection. Section [IV] 
provides a lower bound of not detecting targets achievable 
with these strategies, which is illustrated with a specific 
example in Section [V] Conclusions including future work 
are discussed in Section [VI] Supplementary lemmata used in 
Sections [Til] and [IV] are found in the Appendix. 

II. PROBLEM STATEMENT 

The dependence of detection process on angle occurs 
often in search and detection operations. In general, the 
effectiveness of such an operation also depends on the 
distance between the sensor and the target. However, here, 
the probability of detection is assumed constant as a function 
of range and, hence, the focus is only on the angular 
dependence. For more details on the range dependence, refer 
to Ref. [6]Q 

1 Note that the probability of detection as a function of range is primarily 
a characteristic of the sensor, while the probability of detection as a function 
of angle is primarily a characteristic of the target. 



Fig. 2. Symmetries of the target: Reflection through the short axis of 
the target (dashed line) and reflection through the center of the target (dot- 
dashed line). 



As shown in Fig. Q] the problem is modeled on a two- 
dimensional plan and the observation angle, x, is defined 
as the counter-clockwise angle measured in radian between 
the sensor beam and the short axis of a rectangular (positive 
horizontal axis) target. An observation angle of zero degree 
corresponds to the observation of the long side of the target, 
while an observation angle of tt/2 degrees corresponds to the 
observation of the short side of the target. Targets considered 
will have approximate rectangular symmetries as shown in 
Fig. |2 That is, they possess a reflection axis through their 
short axis (left-right mirror symmetry) and a point reflection 
through their centerQ 

In what follows, the probability of no detection rather 
than the probability of detection is considered; one being 
the complement of the other. Define the single probability 
of no detection as the probability of not detecting the target 
at angle x and denote this single-value real function as g(x). 
Note that the single probability of no detection is even due to 
the reflection symmetry through the short axis of the target 
and periodic due to the reflection through the center of the 
target. Specifically, 

00) = 9{-x), 
g(x) = g(x + ir). 

Next, define the multiple probability of no detection as 
the probability of not detecting a target after n observations. 
Let x be the orientation of the target. Let Ui be the i- 
th angle at which the target is observed relative to x and 
jl = (fio, ■ ■ ■ , Mn-i) be the vector of the n relative ob- 
servation angles. Assume the multiple observation detection 
process is a Bernoulli process, i.e., all observations are 
independent. Then, the multi-observation probability of no 
detection is modeled as the product of single probabilities 

2 Note that the composition of a reflection through the short axis followed 
by a reflection through the center of the target is equivalent to a reflection 
through the long axis of the target (forward/backward mirror symmetry). 



of no detection. In general, however, the exact value of x, 
i.e., the orientation of the target is unknown. To circumvent 
this problem, assume that the target's orientation is uniformly 
distributed and evaluate the average multiple probability of 
no detection by G(jl). Then, 



G(fi) = - 

7T 



dx Yl g( x + [M)- 



(i) 



i=0 



Therefore, the problem amounts to finding search strate- 
gies that minimize G(jT). For simplicity, the probability of no 
detection is taken to mean the average multiple probability 
of no detection in what follows. 

III. A SET OF OPTIMAL SEARCH STRATEGIES 

In this section, a condition ensuring that all partial deriva- 
tives of the probability of no detection are equal to zero 
is derived. From this condition, a set of optimal search 
strategies is then identified and the probability of no detection 
is recast into a form used in the subsequent section. 

Let us first introduce some useful notation and definitions. 
Let i e {0, . . . , n - 1}, Afi = {0, . . . , n - 1} \{i} and j e 
Mi. Define di = d/di. Denote jl* as an optimal point of 
G (jl) . Let a be an integer and b be a positive integer. Define 
the modulo operation as 

a mod b = a - 
Let m be a non-negative integer and /1 = ir/n. Define 

1 n 'f,- 1 

G (m, n) — —l dx j | g(x + mi/i). 



i=0 



Then the following holds. 

Lemma 3.1: The partial derivative can be written as 

diG ($) = — j dx g' (x) j J| g (x + (ij - ^) 



je Mi 



je Mi 

Proof: Apply the partial derivative to the expression 
for G (jl) given by (Q~|). Then 

diG(p)= - I dxg'(x + Hj) J| g(x + Hj). 

n je M z 

Let x — > x— Hi and note that Lemma lATI applies and dictates 
that the integral is invariant under this change of variable. 
Thus, 



diGQS) = 



1 



dx g' (x) Yl 9(x + Hj - Mi)- ( 2 ) 

je Mi 



Because g(x) is even, 



diG (/I) = - 

7T 



dx g' (x) Y[ g(- 

je Mi 



X - fj,j + ^li). 



Let x — > —x and remark that g'(x) is odd. Then 

diG (p,) = -- J dx g' (x) J| g(x - fij + m). (3) 

And the result follows from the average of (|2) and (|3). ■ 

A condition for optimizing G (fl) is thus for the integrands 
of all partial derivatives to be equal to zero. The next Lemma 
identifies a set of search strategies for which this condition 
holds, i.e., optimizes the probably of no detection. 

Lemma 3.2: Define (m, n) to be the search strategy such 
that the separation between two consecutive observations is 
a constant and equal to mir/n with m a positive integer. 
Then the search strategy (m, n) is an optimum of G (fl) and 
defines a subset of all possible optimal search strategies. 
Proof: Evaluate the product in (O at point fl* giving 

n 9{x+n*j — 11*). 

Assume that (m, n) is an optimal search strategy. Then, the 
definition of the strategy (m, n) implies that 

fi* - fi* = m [- (2i - j) + i] fi. 

From definition of the modulo note that 

2i - j 



2i - 



3 



+ (2i — j) mod n. 



Define er,; (j) = (2i — j) mod n and recall that g(x) is 
periodical. Then, 

g(x + fi* — fi*) = g (x - moi (j) fi + mifi) . 

Lemma IA.2I implies that the map (j) is a bijection from 
the set of j to itself. Therefore, 

II 9 {x + fi* - fi*) = \\ g (x - mai (j) fi + mifi) 

je Mt j'e Mi 

= jQ g (x — mjfi + mifi) 

je Mi 

= n -9 ( x - $ + m*) ■ 

je Mi 

And Lemma IXTl then entails that d l G {fl*) — for all i. ■ 

For an optimal search strategy (m, n), it will be useful in 
what follows to cast the integral in the following form. 

Lemma 3.3: Let the search strategy be (m, n). Then the 
probability of no detection is equal to G(m,n). 

Proof: First, remark that only the difference between 
two consecutive observations is specified in Lemma 13.21 
Thus, liberty exists in the choice of the absolute reference 
angle. Without loss of generality, take fiQ to be this absolute 
reference angle. Then, 

fi* =fi* Q + mifi. 

Substituting this expression into (Q]i yields 

1 Tl — l 

G(fi*) = — / dx Y\_9( x + Ho + mifi). 



i=0 



Next, let x —> x — fi*,. Then Lemma I A. 1 1 entails that the 
domain of integration is invariant under such a shift of the 
integration variable. ■ 



IV. A LOWER BOUND OF THE 
PROBABILITY OF NO DETECTION 

In the previous section, a set of optimal search strategies 
was identified, the (m, n) search strategies. In this section, 
a lower bound of the probability of no detection achievable 
with these strategies is proven. 

Let us first introduce some useful notation and definitions. 
Let gcd (•, •) be the greatest common divider. Let r, q and p 
be strictly positive integers such that m — pq, n — pr and 
p = gcd(m, n). Let i s {0, . . . , n — 1}, j e {0, . . . , r — 1} 
and k € {1, . . . ,p}. Define 

h (x) = g (x) g (x + m/i) . . . g (x + (r — 1) mfi) . 



Lemma 4.1: The following identities hold: 

p n—l 

l[h(x + (k-l)fi)=l[g(x + in), (4) 



k=l 



i=0 



g(x + mifi) — h(x) p . 



(5) 



i=0 



Proof: Consider the first identity. Using the definition 
of h(x), write the left-hand side of © as 

p i — l 

n n 9 ( x + + ( k - !) m) • 

k=l j=Q 

The definitions of m and n and of the modulo operation 
imply that 



mj = n 



■ p (qj mod r) 



Using this expression for mj and the periodicity of g (x), 
the right-hand side further becomes equal to 

p i — l 

II J].9(x+p(qj mod r)fi + (k - l)fi). 
k=i j=o 

Next, from Lemma lA~3l the map a (j) = qj mod r is known 
to be a bijection from the set of j to itself. Therefore, the 
product can also be written as 

p i — l 

Y[Y[g(x+ P jfi + (k-l)fi). 

fc=l J=0 

Finally, from Lemma IA.4I the set of (j, k) pairs can be 
mapped to the set of i using the bijection a(k — l,j) — 
(k — 1) + pj. Therefore, 

p r—l n—l 

n n 9 ( x + + ( fc ^ i ) m) = n 9 ^ + ^ ■ 

k=l j=0 i=0 

Now, consider the second identity. From Lemma IA.4I the 
set of i can be mapped to the set of (j, k) pairs using the 
bijection a(j, k — 1) = j + r(k — 1). Therefore, the left-hand 
side of (0 becomes 

p i — l 

IT TT 9 ( x "W* + m (k — l) rpb) . 

k=l j=0 



Next, note that the definitions of m and n imply rar = nq 
and that the definition of /i implies n/x = tt, from which 
follows that m(k — l)r/i — (k — l)g7r. This equality and 
the periodicity of g (x), then imply that 

g (x + mj/i + m(k — 1) rfi) = g (x + rnjn) ■ 

Finally, from the definition of h (x), 

p r—l p 

n n 9 (x + mjn) = n = h( X r. 

k=l j=0 fc=l 

■ 

Theorem 4.2 (Lower Bound Estimate): For any (m, n) 
search strategies, 

G(m,n) > G(l,n). 
Proof: Consider a given (to, n) pair. Then, Lemma |4~T1 
allows to write 



G(l,n) 



i n 



TT 

1 /" = 



dx Y[h(x+(k-l) /it), 



fe=i 



G(m,n) = — dx h(x) . 

TT I x 

Next, let that 2 6 {0, . . . ,p — 1} and define 

p-l 

Xl(x) = Y[h(x + (k-l) f i). 



fc=i 



Assume that 



dx Xi(x)h(x) < / dx Xi+i (x) h(x) 



i+i 



Then 

G(l,n) = - I 2 dx A (x) < - / </./■ A L (.r) //(..") .. 



< 



< 



7T 

1 

7T 

1 /" 7 

7T 



c?x A; (x) h(x) < 



dx A p ~i (x) h(x) 



" 1 = i r dx /i(xf 



= G (to, n) . 



To show that the assumption holds true, proceed by 
generating a partially ordered set. Let Si = (p — I — 1) \i and 
note that \ (x) = A; +1 (x) h(x + 5). Then, the first element 
of the partially ordered set is 

dx A; + i (x) h(x + S)h(x) 1 . 

And the last element is 



l+i 



dx A/+i (x) h(x) 



For greater clarity, the indice of A/ + i (x) and Si are 
suppressed in what follows. Next, apply recursively t 
times h(x) a h(x + S) b < ^h(x) a ^ b h(x) 2a + h(x + 5)' 



from the first element. After summing the resulting geometric 
series, the element generated is 



dx A (x) 



1 



h(x) 



i+i 



(1) h(x) l+1 ~ 2t h(x + 5f 



Remark however that beyond a recursion step defined such 
that 2 i+1 > I + 1 > 2*, the power of h(x) becomes negative 
in the second term. This is problematic because the last 
element of the partially ordered set has only positive powers 
of h(x). To circumvent this problem, let x — > x — 5 in 
the the second term of the element. Since the domain of 
integration remains unchanged (Lemma IA.lt . h(x) is even 
(Lemma IA.5t . and A(x) = A(x — 5) (Corollary IA.6t . then 
the second term is also equal to 



dx A (-x) h(-x) l+1 ~ 2 ' h(-x + S) 



Finally, letting x — > — x, the t-element becomes 
ii 'ii*" 

dx A (x) 



h(x) 



i+i 



l \) h(x + S) l+1 - 2t h(xf 



And the power of h(x) is now greater than or equal to that 
of h (x + S) in the second term allowing again the recursive 
application of the inequality. 

The algorithm thus loops through successive recursion 
steps and changes of variable. Note that for 2 + 1 = 2', 
the last element is generated after the first iteration. Lemma 
IA.7I states that this is the only case for which the algorithm 
terminates in a finite number of iterations. For all other cases, 
the last element arises by letting the number of iterations go 
to infinity. ■ 

V. AN ANALYTICAL EXAMPLE 

In the previous section, a lower bound of the probability of 
no detection was found for an arbitrary g(x). In this section, 
the probability of no detection is evaluated and the inequality 
of Theorem 14.21 is explicitly verified when g (x) = sin (x) 



Lemma 5.1: Let g (x) = sin (x) . Then 

2P(2p- 1)!! 



G (to, n) ■■ 
Proof: Recall that 

i n 



4«p! 



(6) 



p i — l 

G (to, n) = — l dx .9 (x + jm/i). 

77 "'""f k=l j=0 

Then Lemma IA. 81 implies that 

G(m,n) = - / ^dx Y[Y[g[x + — 

J -f k=lj=0 v 



Let g (x) — sin (x) 2 and use the identity [15]: 



sin (rx) = T 1 ]^[ sin I x + — 

3=0 



Then, 



If* 1 
G(m,n) = - I dx J| — T sin (rx) z 



2" fc=l 



1 /" = 



1 



— / da; — -, — rvsin (rx) 
7T /_„ 4p(» - -i) v ; 



2p 



Finally, carrying out the integral [16] and recalling that n 

pr, 

~ 1 (2p-l)!! 

Since 2p!! = 2 p p!, © follows. 

Corollary 5.2: For p = gcd (1, n) = 1, 

2 



G l,n) 



4« 



And Theorem 14.21 follows since (2p — 1)!! = 1 x 3 X • • • X 
(2p- 1) > 1 x 2 x ••• xp=p\. 

VI. CONCLUSIONS AND FUTURE WORK 

In this paper, the angular dependence of the detection 
process which is often overlooked for search and detection 
mission is explicitly accounted for by assuming that the tar- 
get possesses rectangular symmetry. One major consequence 
of this approximate symmetry is that the long side of a target 
is endowed with the largest cross section, which results in 
the highest probability of detection given that the target is 
observed only once. However, since the orientation of the 
target is in general unknown, there is likelihood that it will 
be imaged on the short side, i.e., the smallest cross-section. 
Therefore, the probability of not detecting the target may not 
be zero even if the search area is entirely covered. Making 
several observations of the target in order to increase the 
change of observing its long side is one way to address this 
problem. 

Assuming that the observations are independent, an opti- 
mal search strategy is then to observe the target such that 
the separation between two consecutive observations is a 
constant and equal to a multiple of 180 degrees divided by 
the number of observations. The resulting tactic is simple, 
intuitive and robust (as no prior knowledge of the target 
orientation is required). For example, two observations sep- 
arated by 90 degrees or three observations separated by 60 
degrees will minimize the probability of no detection. 

Having shown that one of these search strategies leads 
to a lower bound of the probability of no detection, work 
is currently underway to prove that it is also a globally 
minimal search strategy. Another logical extension of this 
work is to relax the assumption that information provided 
by subsequent observations is uncorrected (i.e., follows 
a Bernoulli process) as it entails that the probability of 
not detecting a target decreases with increasing number of 
observations even if the observations are co-linear. This is 
questionable as no additional information is gained. 
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Appendix 
Supplementary Lemmata 

Lemma A. 1 (Rotational Invariance): Let w £ R and 
h(x) = h(x + 7r). Then 

dx h(x — uS) = / dx h(x). 

Proof: Let x —> x + uj on the RHS and break the 
integration interval of the resulting integral into [— § — — § [ 
and [— ■§, § — (J\. Then let x — > x — ir in the integral over the 
first interval and recall that by assumption h{x) is periodic. 

■ 

Lemma A.2: Let i G {0, ...,n — 1}, Mi = 
{0,...,n- and j G M l . Then a, (j) = 

(2i — j) mod n is a bijection from Mi to itself and 
its own inverse. 

Proof: Composition gives 

°~i fa (j)) = (2i — (2i — j) mod n) mod n. 
Use the definition of the modulo twice to give 

(2'-i) 



(2i-i) 



3 + 



Recall that j G Mi and note that 



< 



Then 



(2< 



n \ = [^\ n and ^ fa C?')) = 3. ■ 

Lemma A.3: Let r, q be positive integers such that 
gcd (r, q) = 1, i.e., r and q are co-primes. Let i G 
{0, . . . , i — 1} . Then the map er(i) = qi mod r is a 
bijection of the set of i to itself. 

Proof: Proceed with a proof by contradiction. Assume 
this map is not a bijection. Then there exists a pair u, v G 
{0, . . . , r — 1} such that u ^ v and qu mod r — qv mod r. 
Next, assume that u > v then q(u — v) mod r = qw mod 
r = where w = u — v. This implies that qw = ry with y a 
positive integer, as w, q > 0. Because q and r are co-primes, 
i.e., gcd (r, q) = 1 then w = rz with z > 0, which leads to 
a contraction us w = u — v < r — 1. ■ 
Lemma A.4: Let it G {0, ... a — 1}, v G {0, . . . ,6 — 1}, 
and w G {0, . . . , n — 1} where n — ab. Then a (u, v) — 
u + av and er -1 (w) = (w mod a, I — J ) are bijections. 
Proof: Composition gives 

w 



{a-\w)) = 



mod i 



w, 



where the last equality follows from the definition of the 
modulo operation. Similarly, 



a 1 (a(u,v)) = 



(u + av) mod a, 



u + av 



a 



= (u,v) , 



where the last equality follows from the definition of the 
modulo operation and since u < a. Therefore, a(u, v) and 



a 1 (w) are both one-to-one, onto, and inverse of each other. 

■ 

Lemma A.5: h (x) = h (—x). 

Proof: Because of the periodicity of g (x) and nq = 
mr, g (x + mjfi) = g(x+m(j — r)fi). Let j —j + r. Then 
h(x) = g (x) . . . g (x — (r — 1) jfi) and the proof follows 
since by definition g (x) is even. ■ 
Corollary A.6: Consider A/(x) and 5i = (p — k — 
Then A; (x — S{) = h (x — {p — I — 1) fi) . . . h (x) and 
A; [x — Si) = A; (— x) since h(x) is even by Lemma |A. 51 

Lemma A.7: Let i and c be positive integers. Let {^} be 
the set of non-negative integers such that 2 ti+1 > c > 2*\ 
And let {ui} be the set of non-negative integer such that 
Ui = c — 2 4i Ui_i and uq — 1. Then, a fix point exists if and 
only if c = 2 h . 

Proof: After i iterations, 

Ui = c - 2'*i i i_i = c-2*< (c - 2 t - 1 u. i _ 2 ) 
= (c-2* (...(c-2 1 )...)) 
= c [1 - 2 4i (. . . (l-2 t2 ) ...)]+ (_i) i 2*'+-+*i. 

Then ttj = implies that the prime factorization of (fc + 1) 
must be 2° where a is a non-negative integer. Because 
2 tl+1 > (Z + l) = 2 a > 2* 1 , then a = t a and Uj = 
for i > 0. ■ 



Lemma A. 8: 

r- 1 

J=0 



mjTT 



n 



= 3 1 + — 



J7T 

r 



Proof: Recall that r, g and p are positive integers such 
that to = pg, n = pr, and p = gcd(m, n). Then 



9 



TOJ7T 



Use the periodicity of g (a;) and the definition of the modulo 
to further write the right-hand side of the equality as 



g \x + 



r 



r 



TT^j = g [x + (qj mod r) —^j 



Lemma lA3l implies that map a (j) = qj mod r is a bijection 
from the set of j to itself. ■ 
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